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Abstract 

In this paper the dependence of meson-nucleon-nucleon vertex form factors 
is studied as a function of termperature. The results are obtained starting 
from a zero temperature Bonn potential. The temperature dependence of the 
vertex form factors and radii is studied in the thermofield dynamics, a real- 
time operator formalism of finite temperature field theory. It is anticipated 
that these results will have an impact on the study of relativistic heavy-ion 
collisions as the critical temperature for the phase transition from hadronic 
to quark-gluon system is approached. 
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I. INTRODUCTION 



In-medium properties of hadrons and their interactions is a field of high current inter- 
est. For studying phase transitions and thermodynamical properties of the nuclear system 
under extreme conditions, it is essential to determine the temperature (T) and density (p) 
dependence of the nuclear force. Below a critical temperature T c , meson exchange picture 
of nuclear physics provides the natural description for the nucleon - nucleon (NN) interac- 
tion. Quarks and gluons are certainly present in all NN interactions, but it is not always 
necessary to take them explicitly into account especially below a critical temperature T c , 
when the transition into a quark gluon phase takes place. At distances smaller than a fermi, 
the inner structure of hadrons is probed since it involves the short distance or the high mo- 
mentum components of the wave function. This structure is typically taken into account by 
including vertex form factors. In hot and dense matter, the structure of hadrons undergoes 
changes which should lead to a modification of the meson - nucleon coupling constants or, in 
general, the form factors. The temperature [1] and density dependence [2] of some coupling 
constants have been investigated recently, while that of form factors is still fairly unknown 
and will be studied in the present paper for the first time. For this purpose we shall use the 
formalism of Thermo Field Dynamics (TFD) in order to obtain the required temperature 
and density dependence. 

The Thermo Field Dynamics (TFD) which was first suggested by Takahashi and 
Umezawa [3] is a real time operator formalism of finite temperature quantum field theory. 
In this framework, all the operator formalism of quantum field theory at zero temperature 
can be extended directly to finite temperature and density. Therefore in TFD it is possible 
to continue using Wick's theorem and the Feynman diagrammatic approach as in the case 
of zero temperature theory. An auxiliary field is introduced which leads to a field doublet; 
the propagators become 2*2 matrices, the Feynman rules are now algebraic operational 
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rules in the space of 2 * 2 matrices. Recently, some attempts to study in medium properties 
of hadrons at finite temperature have been obtained using TFD [4, 5] . 

A study of relativistic heavy-ion collisions [6] is expected to lead us to a proper equation 
of state of hot and dense nuclear matter. The behaviour of nuclear matter as a function 
of temperature and density is relevant for investigation in nuclear physics, astrophysics, 
cosmology and particle physics. The phase transition from hadronic system to a quark- 
gluon plasma phase and the subsequent hadronisation of the quark-gluon plasma phase can 
provide information on the asymptotic freedom as well as confinement behaviour in Quantum 
Chromodynamics (QCD). Therefore it is important to study theoretically the behaviour of 
the hadronic system at high temperatures. The present report attempts to study such a 
behaviour. In sec. II, a brief review of the TFD is given. In sec. Ill, we calculate the 
renormalization of meson nucleon vertices using the Feynman propagators from TFD. The 
results of calculations and their discussions will be presented in sec. IV. 

II. FINITE TEMPERATURE FORMALISM FOR MESON AND NUCLEON 

PROPAGATORS. 

The Thermo Field Dynamics is a real time operator formalism of quantum field theory 
at finite temperature. The main feature of the TFD is that thermal average of operator A is 
defined as the expectation value with respect to a temperature dependent vacuum, |0 (/?)), 
which is obtained from the regular vacuum by a Bogoliubov transformation. Therefore, we 
have 

< A >= Ti(Ae-^ H -^)/Ti(e-^ H -^) = (0((3)\A\0((3)) (1) 

where f3 = l/ksT with ks being the Boltzmann constant, H is the total Hamiltonian of 
the system, and \i is the chemical potential. Due to the doubling of all degrees of freedom, 
every field is represented in the thermal doublet form in TFD. For example, the thermal 
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doublet for nucleoli field is 



^ a \x) = 



ip(x) 



(2) 



where ip(x) is the ordinary nucleon field and ip is the doublet partner of ip(x). As a conse- 
quence the thermal propagator of a fermion field is a 2 * 2 matrix defined as 



iSp b \ Xl ,x 2 ) =< 0(P)\T[^ a) (^ {b) (^)}\0(P) > 



(a), 



(3) 



The free thermal propagator Sp to be used in perturbation theory can be separated 
into two parts, i.e. the usual Feynman part So, propagator at zero temperature and the 
temperature dependent part St such that Sp b ^ = Sq^ + . For a fermion with mass M, 
the propagators have the form 



S { Q ab) = (^ + M) iab) 



4 ab) = 2mS[f -M 2 )(£> + M) a f : 



f G (p 2 ) ^ 

V G o(P 2 ) J 
sin 2 9 Po \ sin 29 Po ^ 
y \ sin 29 po - sin 2 9 P0 j 



(4) 



(5) 



with 



0(po) 



+ 



O(-Po) 



COSdpo (l + e-)V2 1 (l + e x)i/2' 



(6) 



sin0(p o ) 



e- x ' 2 9{p ) e x / 2 9(-p ) 



(l + e -*)i/2 (l + e*)i/2' 
where a and 6 are Dirac indices and x = f3{po — //). The propagator of the scalar particle 
with mass mg is given by 



A B (p) = A*(p) + A T (p), 



A B 



D (P 2 ) 






-D* (p 2 : 



I 



(7) 



(8) 



with 



D T = -2iri5(p 2 - m%) 



■ Po |sinh20 po 



V 2 



sinh 2 

|sinh20 po sinh 2 0. 



(9) 



Po j 



cosh</> 



1 



sinh0r 



'l- e -\y\)i/2- 
e -|»l/2 



(10) 



Po (1 - e -lvl)i/2' 

where y = (3po. Note that the propagator of vector mesons would be similar to that of the 
scalar mesons except for an additional factor of {—g^v + k fJi k u /m 2 ). 

III. VERTEX FORM FACTORS 



The OBE model of the N-N interaction [7] includes exchange of several mesons. Here we 
consider the meson- nucleon interactions of n.,a.,u- and p-mesons: 



(11) 



c 



pNN 



Here Qbnn are the meson-nucleon coupling constants and f W NN and f P NN are the tensor 
coupling constants for oj- and p- mesons respectively. 

The corresponding meson nucleon form factors are usually defined as: 
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(N(p')\T^\N(p)) = G nNN (t)u(p')i l5 T a u(p) 



(N(p')\UN(p)) 



G aNN (t)u(p')u(p) 



(N(p')\rZ\N(p)) = u{p') <y"G uNN (t) + 



(12) 



(N(p')\T^\N(p))=u(p')r a ^G pNN (t) + ^-F pNN {t)a^q u u{p) 



where G nNN (t), G aNN (t), G uNN (t) (F uNN (t)) and G pNN (t) (F pNN (t)) are form factors of 
7T— , a—, lu — , p— iViV vertices respectively; t = q 2 = q 2 — <f = (p — p') 2 is the 4-momentum 
transfer and M is the nucleon mass. 

In order to investigate the T dependence we calculate three line vertex correction, as is 
illustrated in Fig. 1 for e.g. G n NN(t,T). In accordance with Feynman rules this may be 
written as: 



where A,B- denote 7i,a,u,p mesons, e.g. = G^NNit), a = p' — k,b = p — k. Here 
Sf(o) and A B (k 2 ) are the in - medium nucleon and meson propagators respectively. In 
the framework of TFD the thermal propagator has a 2 * 2 matrix structure, but only the 
11-component refers to the physical field. So, we may use the following representation: 



r A (t,T) = r A (t) + ^A AB (t,T), 



B 



(13) 
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S$(a) = S F (a) = (a + M)[G (a) + G T (a)} , 
A 11 = A(k 2 ) = D (k 2 ) + D T (k 2 ) , 

(14) 

D (k 2 ) = , D T (k 2 ) = -2mn B (k)5(k 2 - m 2 ) , 

k z — m z + ie 

G (a 2 ) = 2 * - . , G T (a 2 ) = 2m5(a 2 - M 2 )N F (a) , 
a — M -\- is 

where N F {a) = 9(a )n F (a) + 0(— a )n F (a), 6>(a ) is the step function, m is the mass of 
corresponding meson and 

UF ^ = efKaoL) + 1 ' fb F (ci) = e/3(|ao| ^ + 1 , n B {k) = -^- - (15) 

are the fermion, antifermion and boson distribution functions respectively. It is clear that 
when one substitutes Eq. (14) into Eq. (13), A AB (t, T) will be separated into two parts: 
one refers to the naive zero temperature contribution, and the other one depends on T and 
p [4]. We shall concentrate on the latter part, which may be rewritten as follows: 

r rlk 4 

A AB (t,T) = J —T B (k 2 )(a + M)T A (t)(b + M)T B (k 2 )x 

(16) 

x [G (a)G (b)D$(k 2 ) + 2G (a)G T (b)D*(k 2 )} 
Note that we have neglected terms that are quadratic in the temperature dependent distri- 
bution function. Now, using (16) and (14) in (13) we get the following expressions for the 
form factors: 

-—W AB (t, k 2 ,T) 

(27T) 1 

x[G (a)G (b)D^(k 2 ) + 2G (a)G T (b)D*(k 2 )] 

where the explicit formulas for WA B {t, k 2 , T) may be found in the Appendix. Here we write 
expressions for W un (t, k 2 , T) and W pn (t, k 2 , T): 



dk^ 

u W)*''^"^' (17) 
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w^(t,e,T) = 3 ^lmm j 2[2M * - (a6 )] - ; 

(18) 

W p „(t, k\T) = Z^mm | 2[2M2 _ {ab)] _ ; 

where W^{t, k 2 , T) and W^t, k 2 , T) denote the contribution from the ir- exchange diagram 
to the G^NNit.T) and G P NN(t,T) respectively. 

IV. RESULTS AND DISCUSSIONS 

In order to start the calculations, a set of free space meson - nucleon form factors are 
chosen. We choose the OBE monopole form factors [8] (Bonn A): GsNNit) = Qbnn{.^b ~ 
m 2 B )/(A 2 B — t), where , for example, g 2 NN /^ = 14.09 , and A nNN = 1005MeK. Let us first 
discuss the T dependence of meson - nucleon coupling constants - Qbnn{T^) = GsNNit = 
171%, T). The variation of the coupling constants g n NN(T) / g n NN(T = 0), g a NN(T) / g a NN(T = 
0), guNN{T)/gu,NN{T = 0) and g pN N(T)/g pNN (T = 0) with temperature are displayed in Fig. 
2. It is clear that they are nearly independent of the temperature below T C B , then change 
rapidly for T > T C B . A similar behavior was also predicted by Zhang et al. [4] and by 
Dominguez et al. [1] for pion - nucleon coupling constant g^NNiT)/ g nNN (T = 0). But here 
in Fig. 2, an unexpected result is that, the oj — N coupling constant, guiNNiT) increases while 
the coupling constants of all other mesons decrease! Let's consider a possible origin of this 
controversy, by comparing the medium modification of uNN and pNN coupling constants. 
In the present model the modifications arise from the triangle diagrams (Fig. 1). Actual 
calculations show that, for any GBNN(t,T) the triangle diagram with pion exchange gives a 
dominant contribution. It is clear from Eq. (18) that at small t the explicit expressions for 
gu)NN{T) and g pNN (T) formally coincide. However, W pw = — ^f L - and hence the two have 
an opposite sign. The factor (—1/3) arises from the isotopic spin structure. In fact, using 
Eq. (13) and Fig.l, g^NN^T) and g P NN{T) may be written in a schematic way: 



9u>nn(T) guNN{T — 0) + g^jvjv X! tp9u>nnT/3 + ■ ■ ■ 

P 

(19) 

g P NN{T)r a m g pN N{T = 0)r a + gl NN ^ TpT a g pNN Tj3 + ... 

where a and (5 are isospin indices. We omit the spin variables, since uj and p have the same 
spin structure. Now, using = 3 and T^T a Tp = —r a , it is clear that the contribution 

of the leading triangle diagram with pion exchange to g U NN(T) and g P NN(T) has different 
signs. In other words, the term A pn in Eq. (13) is negative, while is positive. 

In general, we conclude from the results (Fig. 2) that at a certain critical temperature T c s , 
the couplings gsNN change drastically. The thermal behavior of g^NN^T) = G n NN(t,T)\ t ^o 
has been investigated earlier [1] and we get here a similar behaviour of g^NN^T) as a function 
of temperature. The temperature T c , where g n NN(T) changes dramatically, was interpreted 
in [1] as a signal for the quark - gluon deconfinement phase transition. Indeed, near T c 
the associated mean square radius < r 2 BNN >= 6[J^ In Gbjvat(^, T)]| t _to is a monotonically 
increasing function of T and in fact diverges at the critical temperature. A similar behavior 
of < 

r BNN > 1S a l so found in our calculations and is illustrated in Fig. 3. As the critical 
temperature, T c s , is approached the strength of coupling of 7r, a and p mesons to nucleons is 
quenched, at the same time, the size of nucleons as probed by the appropriate meson gets 
bigger. However, it is to be stressed that T c s is not the same for all mesons: T£ ^ 360MeV, 
T° ss 95MeV, ^ 175MeV and T£ ^ 200MeV. In OBE picture this means that , in 
the temperature region e.g. 200MeV < T < 300MeV, the p meson exchange is no longer 
important while the pion exchange is still important. On the other hand, the a and n mesons 
are mainly responsible for the attraction between nucleons. So, the quenching of the aNN 
coupling constant at T > T c °" leads to the vanishing of the bound state, as it was predicted 
earlier [4]. 

The density dependence of vertex form factors has been studied in greater detail [2, 8]. 

9 



It was found that most of the vertex form factors are quenched at high densities and it was 
anticipated that the temperature dependence is likely to yield results that are qualitatively 
similar to those of density dependence. Now, let's consider their temperature dependence in 
some detail. The form factors as a function of momentum transfer at several temperatures 
are displayed in Fig. 4. It is seen that, the in- medium effects lead to the suppression of 
GwNN(t,T), G(rNN(t,T) and G P NN(t,T). The temperature dependence of GuNN(t,T), is 
opposite to that of G p nn due to it's isotopic structure as outlined above. 

For practical calculations a parametrization of these GsNN(t,T, p) form factors is needed. 
At small momentum transfer we can parametrize them by a monopole form : 

G BNN (t,T,p) = g B (T,p)(A 2 B (T,p)-ml)/(Al(T,p)+t) (20) 

where in general the effective mass of a meson, tub, is also temperature dependent. Consid- 
eration of this dependence is beyond the scope of the present paper. It will be carried out in 
a self consistent way, taking into account both the self energy and the triangle diagrams, in a 
forthcoming paper. But here, for simplicity, one may consider this as just a parametrization 
and choose the parameters g B (T,p) and A B (T,p). Their temperature dependence is still 
unknown. Here the T dependence maybe represented, for T << T c s , in a polynomial form 



as: 



" b{ LP) 1 - a B ( P )(T/T c B f - (3Up)(T/T c B r 



g B (T = 0,p = 0) 

(21) 



l-a B (p)(T/T c »y-f3 B (p)(T/T c »y 



A B (T = 0,p = 0) 

The calculated form factors are fitted to this form and the paramters a and (3 are determined. 
The results are presented in tables I, II, and III for densities p = 0, p = po and p = 3po, 
respectively (with p = 0.17/m -3 - the density of normal nuclear matter). The results show 
that the density dependence of parameters, and hence, the form factors is rather weak. 
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What is the reason for such a result? The answer may be found by considering expression 
in Eqs. (14) - (17). It is clear that the density, p, influences only through the chemical 
potential \i defined from the particle and antiparticle distribution functions as 



where p is involved in Eq. (15) in the exponent. The argument of this exponent is /3(|a |±p). 
Since |a | ~ M = 939 MeV is much larger at moderate densities than p (p AOMeV for 
p = 3po) T = 150MeV), the result is weakly sensitive to p. 

In summary, we have considered the temperature dependence of meson - nucleon form 
factors and coupling constants in TFD formalism. It is shown that at a critical temperature, 
where the coupling constant changes drastically and the associated mean square radius 
diverges is different for different mesons. The temperature dependent vertex form factors 
are parametrized in a simple monopole form and the T-dependence of these parameters is 
clarified. These form factors may be used in the calculation of the in - medium NN cross 
sections [9] and in investigations of the properties of hot dense matter. 

In the present calculations, the monopole form factors with the parameters (at T = 0) 
given by the Bonn group, [7] are used. But the next question which may arise is: are the 
results sensitive to the shape of the input form factors? To answer this question a different 
set of form factors obtained by Meissner et. al., [10] in the framework of a topological 
soliton model have been used. The thermal behaviour of such form factors is shown in Fig. 
5. Now, comparing Fig. 4 and Fig. 5 it is clear that only the aNN form factor is affected 
the most. The reason is that even in free space and at T = the aNN form factor in 
the topological soliton model [10] is much harder (A CT ~ 2700MeV) than those of the Bonn 
potential (A CT ss 2000MeV). Nevertheless our main conclusions about thermal behavior of 
different meson nucleon vertices remain valid. 

The variation of the meson-nucleon coupling constant with temperature (Fig. 1) clearly 




(22) 
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indicates that the attractive NN force due to a-exchange decreases quite rapidly while the 
repulsive N-N force due to u;-exchanges increases. This would lead to the fact that the 
nuclear matter is quite likely un-bound at high temperature. In fact it will look more like 
a hard sphere gas since the repulsive interaction at short distances will dominate. The 
attractive interactions in such a case play only a small perturbative role. The hot and dense 
nuclear matter could be approximated as a free gas with an excluded volume around each 
nucleon. Perhaps the extensive studies [11] along this line can be justified on the basis of 
the results obtained in this paper. The nuclear matter would be hard to compress. 

It would be quite interesting to know, if the critical temperature, Tf , maybe uniquely 
considered as a point of transition from the hadron state to a quark gluon state, as it was 
suggested [1]. The usual way of estimating T c is based on QCD lattice calculations [12] 
or temperature dependent quark - gluon potentials [13]. In such calculations it is usual 
to start from high temperatures and then decreases the temperature, looking for a phase 
transition from the quark gluon state to a hadron state i.e. hadronization of the quarks. 
Clearly this method gives a unique T c , since at any stage there is either a hadron or a 
quark gluon plasma phase. Contrary to such studies the present calculations explore the 
transition froma hadron to a quark-gluon phase and find that the critical temperature for 
hadron — > quark-gluon transition depends on the kind of hadrons under consideration. A 
similar result, where hadrons and quark - gluon plasma coexist have been obtained in lattice 
calculations [12]. In conclusion it is important to emphasize that the dynamics of hadrons 
and of quark-gluon plama at finite temperature is poorly understood. The question of critical 
phenomenon, in particular critical temperature, in the hadronic systems is not well-known. 
Studies in nuclear matter are needed urgently to clarify the critical transition from a hadronic 
system to the quark gluon phase which will depend on both density and temperature. The 
experiments with relativistic heavy ion colliders would need such an understanding to clarify 
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the production of quark-gluon plasma from hadrons and the eventual hadronisation of the 
quark-gluon system. It is possible that there will be clear and unambiguous signals for the 
formation and subsequent hadronisation of the quark-gluon plasma [14]. 
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APPENDIX 



Here explicit expressions for Wab^, k 2 , T) introduced in eq. (7) are given in detail. 



W vn = G n {k 2 )[M 2 - {ah)} ■ W na = W^\ G ^ Ga ■ 



W nu> = A^\M 2 - (ab)} + 3F?(k 2 )F?(k 2 )z ; W vp = -W v 



K = [Fr(fc2)]2[4 _ j_] + 3^ [FM )] , 



m 



7T 7T 



-3<? 2 (fc 2 )[M 2 + («&)], 



g W<77r I G,r— >G CT ) 



W w = [M 2 + (a&)]A£ ; W,, = 3W^U^p ; 
K = -\F?{k 2 )] 2 [A - 4] + ^ 2 ^^ 2 ^ • 

^ = ^^{2[2M 2 -(a6)] 



4M 2 
3i^(*)t" 



1 



2Ff(t) 



W wu = -W^\G^G a ; 



- 2T^[2M 2 - (ab)} + ^\k 2 \M 2 - (ab)} + 2(ak)(bk)} + 



M 2 



2M 2 Ff{t) 



[3rft + r%(k 2 t- z 2 )}+ 



+3zF% - 

+ ^l z K ab ) + 3M '] - 2(ak)(bq) - 2{aq){bk)\ \ ; 



r - = [Fr(A; 2 )] 2 [2-^ 2 f ] + — [F ^ (A;2)]2 
2M 2 

[^ 2 )] 2 -ra 2 ; 



4M 2 



1 2 — 



r- = Fr(A; 2 )F-(A; 2 ) , 



w, 



per 



-w, 



pir | Gtt—>G& i 



w, 



g M^a;7r|aj— >7r 5 



pec" 

i 



pp 



-w 



LOU) \U)^p ) 



where G n = G nNN (k 2 ), G a = G aNN (k 2 ), F" = G wNN , F£ = F^nn , z = (kg) and q 2 = t. 
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FIGURE CAPTIONS 



Fig. 1. Feynman diagrams for pion - nucleoli vertex, The solid line is for nucleon. Dashed, 
dotted, dot - dashed and wavy lines for n, a, uj and p mesons. 

Fig. 2. The ratio of meson- nucleon coupling constants at finite temperature T and at T = 
as a function of temperature, T. 

Fig. 3. The ratio of mean square radius at finite temperature T and at T = as a function 
of temperature, T. 

Fig. 4. Meson nucleon form factors at several temperatures as a function of q 2 using 
parametrisation of BONN group [7]. a) irNN vertex; b) aNN vertex; c) uNN vertex 
and d) pNN vertex. 

Fig. 5. Meson- nucleon form factors at several temperatures as a function of q 2 using 
parametrisation of Meissnar et al [10]. a) irNN vertex; b) aNN vertex; c) uNN 
vertex and d) pNN vertex. 
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TABLES 



TABLE I. Parameters of vertex form factors in Eq.s (20)-(21) at zero density p = 0. 



Meson 


y 


Pa 


a\ 


Px 


T c 


7T 


-0.201 


0.884 


0.6237 


-0.08419 


w 360 


a 


0.3252 


0.8381 


-0.3245 


0.9184 


» 95 


P 


0.4384 


0.03194 


1.198 


-0.5245 


» 200 


UJ 


-0.8951 


-0.2235 


0.7322 


-0.5258 


« 175 




TABLE II. 


The same as 


in Table I , but for p = 


Po 




Meson 


y 


Pa 


a\ 


P\ 


T c 


7T 


-0.1966 


0.8899 


0.6181 


-0.0750 


w 360 


a 


0.3254 


0.8378 


-0.3301 


0.9329 


w 95 


P 


0.4379 


0.03347 


1.2734 


-0.6071 


« 200 


UJ 


-0.8966 


-0.2168 


0.7332 


-0.5366 


w 175 




TABLE III. 


The same as 


in Table I , but for p = 






Meson 


a g 


Py 


a x 


Px 


T c 


7T 


-0.1868 


0.9012 


0.6056 


-0.05452 


w 360 


a 


0.3258 


0.8366 


-0.3429 


0.9674 


» 95 


P 


0.4379 


0.03347 


1.1934 


-0.5047 


» 200 


UJ 


-0.8994 


-0.2026 


0.7350 


-0.5609 


» 175 
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